Enhancing decoherence times in superconducting qubits via circuit design 
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We study decoherence effects in qubits coupled to environments that exhibit resonant frequen- 
cies in their spectral function. We model the coupling of the qubit to its environment via the 
Caldeira-Leggett formulation of quantum dissipation/coherence, and study the simplest example of 
decoherence effects in circuits with resonances such as a dc SQUID phase qubit in the presence of an 
isolation circuit, which is designed to enhance the coherence time. We emphasize that the spectral 
density of the environment is strongly dependent on the circuit design, and can be engineered to 
produce longer decoherence times. We begin with a general discussion of superconducting qubits 
such as the flux qubit, the Cooper pair box and the phase qubit and show that in these kinds of 
systems appropriate circuit design can greatly modify the spectral density of the environment and 
lead to enhancement of decoherence times. In the particular case of the phase qubit, for instance, we 
show that when the frequency of the qubit is at least two times larger than the resonance frequency 
of the environmental spectral density, the decoherence time of the qubit is a few orders of magnitude 
larger than that of the typical ohmic regime, where the frequency of the qubit is much smaller than 
the resonance frequency of the spectral density. In addition, we demonstrate that the environment 
does not only affect the decoherence time, but also the frequency of the transition itself, which is 
shifted from its environment-free value. Second, we show that when the qubit frequency is nearly the 
same as the resonant frequency of the environmental spectral density, an oscillatory non-Markovian 
decay emerges, as the qubit and its environment self-generate Rabi oscillations of characteristic time 
scales shorter than the decoherence time. 

PACS numbers: 74.50.+r, 85.25.Dq, 03.67.Lx 



I. INTRODUCTION 

The possibility of using quantum mechanics to manip- 
ulate information efficiently has lead, through advances 
in technology, to the plausibility of building a quantum 
computer using two-level systems, also called quantum 
bits or qubits. Several schemes have been proposed as 
attempts to manipulate qubits in atomic, molecular and 
optical physics (AMO), and condensed matter physics 
(CMP), however it is still very difficult to implement a 
scheme that gives both long decoherence times and is 
scalable. 

In AMO, the most promising schemes are trapped ion 
systems^, and ultracold atoms in optical lattices 2 -. On 
the CMP side, the pursuit of solid state qubits has been 
quite promising in spin systems^ and superconducting 
devices^—. While the manipulation of qubits in AMO 
has relied on the existence of qubits in a lattice of ions 
or ultra-cold atoms and the use of lasers, the manipula- 
tion of qubits in CMP has relied on the NMR techniques 
(spin qubits) and the Josephson effect (superconducting 
qubits). Integrating qubits into a full quantum computer 
requires a deeper understanding of decoherence effects in 
a single qubit and how different qubits couple. 

In AMO systems Rabi oscillations in single qubits have 
been observed over time scales of miliscconds since each 
qubit can be made quite isolated from its environment^, 
however it has been very difficult to implement multi- 
qubit states as the coupling between different qubits is 
not yet fully controllable. On the other hand, in super- 



conducting qubits Rabi oscillations have been observed^ 
over shorter time scales (500ns), since these qubits are 
coupled to many environmental degreees of freedom, and 
require very careful circuit design. For superconducting 
qubits, it has been shown experimentally^ that sources 
of decoherence from two-level states within the insulat- 
ing barrier of a Josephson junction can be significantly 
reduced by using better dielectrics and fabricating junc- 
tions of small area (< 10 /mi 2 ). 

In this manuscript, we use the Caldeira-Leggett for- 
mulation of quantum dissipation to analyze decoherence 
effects of generic superconducting qubits coupled to envi- 
ronments that exhibit a resonance in their spectral den- 
sity. This is an extension of our previous work^, where 
the general idea presented in this paper in the context 
of dc- SQUID phase qubits. Here, we show first how 
the characteristic spontaneous emission (relaxation) life- 
times T\ for flu x 10 ' 11 , phased and charge^ qubits can 
be substantially enhanced, when each of them is cou- 
pled to an environment with a resonance, provided that 
the frequency of operation of the qubit is about twice as 
large as the frequency of the environmental resonance. 
Second, we show that the coupling to the environment 
does not only cause decoherence, but also changes the 
qubit frequency. Lastly, we show that when the qubit 
frequency is nearly the same as the resonant frequency 
of the environmental spectral density, an oscillatory non- 
Markovian decay emerges, as the qubit and its environ- 
ment self-generate Rabi oscillations of characteristic time 
scales shorter than the decoherence time. 
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The paper is organized as follows. In Sec. QTJ we de- 
rive the environmental spectral density for flux, phase, 
and charge qubits when each of them is coupled to an 
environment with a resonance. In Sec. IIII| we use the 
Caldeira-Leggett formulation of quantum dissipation to 
derive the Bloch-Rcdficld equations and to calculate the 
decoherence and relaxation times of these qubits. There, 
for the particular case of the phase qubit, we show that 
when the frequency of the qubit is at least two times 
larger than the resonance frequency of the environmen- 
tal spectral density, the decoherence time of the qubit is 
a few orders of magnitude larger than that of the typical 
ohmic regime, where the frequency of the qubit is much 
smaller than the resonance frequency of the spectral den- 
sity. In Sec. IIV| we calculate the renormalization of the 
qubit frequency due to dressing of the two-level system 
with environmental degrees of freedom. In Sec. El we de- 
rive an essentially exact non-perturbativc master equa- 
tion, when these systems are near resonance and where 
the rotating wave approximation can be applied. Our 
results indicate that when the qubit frequency is nearly 
the same as the resonant frequency of the environmen- 
tal spectral density, an oscillatory non-Markovian decay 
emerges, as the qubit and its environment self-generate 
Rabi oscillations of characteristic time scales shorter than 
the decoherence time. Finally, we present our conclusions 
in Sec. ED 



II. SPECTRAL DENSITIES OF 
SUPERCONDUCTING QUBITS WITH 
ENVIRONMENTAL RESONANCES 

In this section, we discuss three examples of supercon- 
ducting qubits coupled to environments (circuits) which 
have resonances in their spectral functions. We begin 
our discussion with a simple flux qubit, then move to a 
charge qubit and finally to a phase qubit. 
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FIG. 1: (Color-online) Flux qubit measured by a dc-SQUID 
gray (blue) line. The qubit corresponds to the inner SQUID 
loop with critical current I c and capacitance Cj for both 
Josephson junctions denoted by the large x symbol. The 
inner SQUID is shunted by a capacitance C 3 , and environ- 
mental resistance R and is biased by a ramping current 
The dc-SQUID loop has junction capatitance Co and critical 
current I c q. 



decoherence and relaxation time scales in such a system it 
is necessary to understand how noise is transferred from 
from the dc-SQUID to the qubit. 

For the circuit displayed in Fig. [TJ the classical equa- 
tion of motion for the dc-SQUID (outer-loop) is 

C o 0+— I cO sin0- — 1+ / dt'Y(t-t')^(t')=0 (1) 

$0 1*0 Jo 

where <j> is the gauge invariant phase across the Josephson 
junction of the dc-SQUID, I c q is the critical current of 
its junction, $o = h/2e is the flux quantum. The last 
term in Eq. ([TJ) is the dissipation term due to effective 
admittance Y(uj) felt by the outer dc-SQUID. In this 
case the total induced current in the dc-SQUID (outer 
loop) isi^ 
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(2) 

Here <p p = (</>i + 02 ) and <f> m — (02 — </>i) are the sum 
and difference of the gauge invariant phases 0i and 02 
across the junctions of the inner SQUID, Ld c is the self- 
inductance of the inner SQUID, and J\ is the bilinear 
coupling between m and p at the potential energy min- 
imum. The term S(j>o = TrM q I c i r /$o, where I C i r is the 
circulating current of the localized states of the qubit 
(described in terms of Pauli matrix a z ), and M q is the 
mutual inductance between the qubit and the outer dc- 
SQUID. 

For the charge qubit coupled to a transmission line 
resonato r 13 ! 15 shown in Fig. [21 the classical equation of 
motion for the charge Q is 



KM 



ui 2 Lj(uj) 



1 
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+ — +iwZ{u) Q(lu) (3) 



Here, V g is the gate voltage, 



C g is the gate capacitance, 
Lj and Cj are the Josephson inductance and capacitance 
respectively, Z(uS) is the effective impedance seen by the 
charge qubit due to a transmission line resonator (cavity), 
and Q is the charge across C g . 
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FIG. 2: Circuit diagram of the Cooper-pair box. The super- 
conducting island (large x) is connected to a large reservoir 
through a Josephson junction with Josephson energy Ej and 
capacitance Cj. The voltage bias V g is provided through a res- 
onator (cavity) having environmental impedance Z{uj), which 
is connected to the gate capacitance C g as shown. 



In Fig. [TJ we show a flux qubit (inner-loop), which 
is measured by a dc-SQUID (outer loop). To study the 



In Fig. [31 we show the circuit for a phase qubit corre- 
sponding to an asymmetric dc-SQUI D 12 i 16 . The circuit 
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elements inside the dashed box form an isolation network 
which serves two purposes: a) it prevents current noise 
from reaching the qubit junction; b) it is used as a mea- 
surement tool. 

The classical equation of motion for the phase qubit 
shown in Fig. [3] is 

0,7+1^0 8*17-^/+ [ dt'Y(t-t')j(t') = 0, (4) 

where Co is the capacitance, I c q is the critical current, 
and 7 is the phase difference across the Josephson junc- 
tion J (large x in Fig. [3]), while / is the bias current, 
and $o = h/2e is the flux quantum. The last term of 
Eq. (fj| can be written as iu)Y(oj)j(u>) in Fourier space. 
The admittance function Y(lu) can be modeled as two 
additive terms Y(uj) = Y iso {uj) + Y int (uj). The first con- 
tribution Yi 80 {ui) is the admittance that results when a 
transmission line of characteristic impedance R is at- 
tached to the isolation junction (here represented by a 
capacitance C and a Josephson inductance L) and an 
isolation inductance L\. Thus, Yi so (u>) = Z~1(oj) where 

Z iso {uj) = (iojLi) + [R- 1 + icuC + (iujLy 1 ]^ 1 is the 
impedance of the isolation network shown in Fig. [3] The 
replacement of the isolation junction by an LC circuit 
is justified because under standard operating conditions 
the external flux $ a varies to cancel the current flowing 
through the isolation junction making it zero biased^. 
Thus, the isolation junction behaves as a harmonic os- 
cillator with inductance L which is chosen to be much 
smaller than L\. The second contribution Yi nt {uj) is an 
internal admittance representing the local environment 
of the qubit junction, such as defects in the oxide bar- 
rier, quasiparticle tunneling, or the substrate, and can be 
modeled by Yi nt (uj) — (Rq + iujLq) -1 , where Rq is the 
resistance and Lq is the inductance of the qubit as shown 
in Fig. O 
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FIG. 3: (Color-online) Schematic drawing of a phase qubit 
with an RLC isolation circuit. The phase qubit is shown 
inside the solid (red) box, the RLC isolation circuit is shown 
inside the dashed box to the left, and the internal admittance 
circuit is shown inside the dashed box to the right. 



written in terms of Pauli matrices Ci (with i = x,y, z) 
and boson operators bk and b\ . The first term in Eq. J5]) 
represents a two-level approximation for the qubit (sys- 
tem) described by states |0) and |1) with energy differ- 
ence ftwoi- The second term corresponds to the isolation 
network (bath) represented by a bath of bosons, where 
bk and b\ are the annihilation and creation operator of 
the fc-th bath mode with frequency LUk- The third term 
is the system-bath (SB) Hamiltonian which corresponds 
to the coupling between the environment and the qubit. 

At the charge degeneracy point for the charge qubit 
(gate charge N g = 0), at the flux degeneracy point (ex- 
ternal flux $ cxt = 7r<f>o) for the flux qubit, and at the 
suitable flux bias condition for the phase qubit (external 
flux $ a = Li4>q) Hsb reduces to 



H SB = \o x K(l\v\Q)^\ hl 
k 



(1 



(0) 



where v = <p f° r the flux qubit, v — Q for the charge 
qubit, and v = 7 for the phase qubit. The spectral den- 
sity of the bath modes 



has dimensions of energy and can be written as 

/$o x 2 
J(w) = LuReY(uj) — 
\ 2tt 

for flux and phase qubits or as 



J(cj) = 2hujReZ(uj) — 

h 



(7) 



(8) 



(9) 



for charge qubits. 

For the flux qubit circuit shown in Fig. [IJ the shunt 
capacitance C s is used to control the environment, while 
the Ohmic resistance of the circuit is modelled by R. In 
this case, the environmental spectral density is^ 



a\u> 



(10) 



when uj m 3> max(f2i, Woi)i an d when the dc-SQUID is 
far away from the switching point to be modelled by an 
ideal inductance Lj. The inner oscillator frequency uj m = 
y/2/ LdcCj and the qubit frequency is cjoi ■ The external 
oscillator frequency 

n x = yj2nif/c a $ [i - (V4 eff ) 2 ] 1/4 , (11) 



The equations of motion described in Eqs. (pj, ([3]), and 
Q, can be all approximately described by the effective 
spin-boson Hamiltonian 



H = 



01 



^ ?iu k blb k + H SB , 



(5) 



and plays the role of the resonant frequency, where I c is 
the critical current for each of two Josephson junctions. 
Also, Ti = \/{C s R) corresponds to the resonance width, 
and 
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reflects the low frequency behavior. The coupling be- 
tween the flux qubit and the outer dc-SQUID emerges 
from the interaction of the persistent current I C i r of the 
qubit and the bias current h of the dc-SQUID via their 
mutual inductance M q . 

The spectral density for the charge qubit shown in 
Fig. [2] is obtained via Eq. ([9]) by determining the real 
part of the impedance Z(u>). In this case, we need to 
solve for the normal modes of the resonator and trans- 
mission lines, including an input impedance R at each 
end of the resonator. This procedure results in the spec- 
tral density^ 



J 2 M = 



r. 



£ C ( w - n 2 y + (r 2 /2) 5 



(13) 



where ^2 is the resonator frequency, I is resonator length, 
c is the capacitance per unit length of the transmission 
line. The quantity Yi = f^/Q where Q is the quality 
factor of the cavity. 

For the phase qubit shown in Fig. [31 the spectral den- 
sity J(w) is given by Eq. (|5J), and can be written in the 
compact form J(u>) = Ji SO (uj) + Ji n t(ui). The spectral 
density of the isolation network is 



Jiso M 



^) (i -w 2 /o 2 ) 2 + 4w 2 r 2 /n 4 ' 



(14) 



where a = L 2 / [(L + LxfR] w (L/Lxf/R is the leading 
order term in the low frequency ohmic regime, 



Q = 




(15) 



is essentially the resonance frequency, and T = 1/ (2CR) 
plays the role of resonance width. Here, we used L± L 
corresponding to the relevant experimental regime. No- 
tice that Ji SO (u}) has Ohmic behavior at low frequencies 
since 



lim 



$o\ 2 / L \ 2 1 



2tt 



LiJ R' 



(16) 



but has a peak at frequency fl with broadening controlled 
by F. In addition, notice that the parameter 
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(17) 



is independent of C. Therefore, when there is no capaci- 
tor (C — > 0), the resonance disappears and 



2tt J l + 4u; 2 r 2 /ft 4 ' 



(18) 



reduces to a Drude term with characteristic frequency 
il 2 /2T ps RJL. The internal spectral density of the qubit 



Jint (w) 



(u/R ) 



2tt J 1 + to 2 Ll/Rl 



(19) 



is a Drude term with characteristic frequency Rq/Lq. No- 
tice that Ji n t{oj) also has Ohmic behavior at low frequen- 
cies since 
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(20) 



Having discussed the spectral functions of the environ- 
ment of three different types of superconducting qubits 
and their corresponding circuits, we move next to the dis- 
cussion of the decoherence properties within the Bloch- 
Rcdfield description of the time evolution of the density 
matrix. 



III. BLOCH-REDFIELD EQUATIONS: 
DECOHERENCE PROPERTIES 



In this section, we investigate the relaxation T\ and de- 
coherence Ti times, for the three types of superconduct- 
ing qubits (flux, charge and phase) described in sectionlHl 
In order to obtain T\ and T 2 , we write the Bloch-Rcdficld 
equation: 



,10 



Pnm — i^nmPnm ~t~ ^ ^ RnmklPkl 
kl 



(21) 



for the density matrix p nm of the spin-boson Hamiltonian 
in Eq. ([5]) and ([6]) derived in the Born-Markov limit. Here 
all indices take the values and 1 corresponding to the 
ground and excited states of the qubit, respectively, while 
^nm = (E n — E m )/h is the frequency difference between 
states n and m. The Redficld rate tensor is 



R 



— _r'- 1 ) — p( 2 ) I X I r p( 2 ) fOOA 

nmkl — 1 Imnt 1 Imnk ' u nk>- T 0/ m l nrrk , {^^) 



where repeated indices indicate summation, and 

dte -i«>n k t ( HsBM (t)H S B,nk (0)) , (23) 



1 Imnk ~ 11 



.(1) 
Imn 

r (2) - h- 2 

1 Imnk ~ 11 
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H S B,lm(0)HsB,nk(t))- (24) 



Under these conditions, the relaxation rate becomes 
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(25) 



where we used the transition probability |(0|^|1)| 2 = 
[Mwoi] 1 to define the variable M, which has dimen- 
sions of mass x area (or energy x squared-time) and is 
refered to as the mass of the qubit. For the flux and 
phase qubits, the mass is M = ($o/27r) 2 Co, where Co is 
the capacitance of the qubit. For the charge qubit, the 
mass is M = He/Io, where Iq is the critical current. The 
frequency woi is the qubit frequency. 

The physical interpretation of T^ 1 is as follows. For 
the system to make a transition it needs to exchange 
energy E = Hujqi with the environment using a single 
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FIG. 4: (Color-online) Ti (in seconds) as a function of qubit 
frequency cjoi . The solid (red) curves describe the phase qubit 
with RLC isolation network (Fig. [3]) with parameters R = 
50 ohms, Li = 3.9nH, L = 2.25pH, C = 2.22pF, and qubit 
parameters Co = 4.44pF, Ro — oo and Lq = 0. The dashed 
curves correspond to an RL isolation network with the same 
parameters, except that C = 0. Main figure (T = 0), inset 
(T = 50mK) with ft = 141 X 10 9 rad/sec. 



boson. The factor coth(hLo i/kRT) = n(cjQi) + l+n(u!Qi) 
captures the sum of the rates for emission (proportional 
to n(woi) + 1) and absorption (proportional to n(cJoi) of 
a boson), where n(idoi) = [cxjp(hujQi/kBT) — 1] is the 
Bose function. 

Similarly, the decoherence rate given by the off- 
diagonal elements of the reduced density matrix p is 
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where the dephasing rate 



- = |(0|H0)-(i|Hi: 



f lim ^l 2 k B T 

w->0 UJ 



(26) 



(27) 



This contribution originates from dephasing processes 
which randomize the phases while keeping energy con- 
stant, i.e., transitions from a state into itself. Hence, 
they exchange zero energy with the environment and 
J(0) enters. The prefactor measures which fraction of 
the total environmental noise leads to fluctuations of the 
energy splitting, such that only the components of noise 
which arc diagonal in the basis of energy cigenstates leads 
to pure dephasingi^. This effect, is thus complemen- 
tary to that of the off- diagonal transition matrix ele- 
ment (0|i/|l) which appears in T\. The zero frequency 
argument is a consequence of the Markov approxima- 
tion. Typically, T^ 3> T\. This can be seen for ex- 
ample, in the Hamiltonian of the phase qubit, where 
(0|7|0) ~ (1|t|1) ~ d 3 U/d'y 3 is the cubic correction to 
the potential of the phase qubit. 

In Fig. [4] Ti is plotted for the phase qubit as a function 
of the qubit frequency woi in the case of spectral densities 
describing an RLC [Eq. (fT3))] or Drude [Eq. (JT5J)] isolation 
network at fixed temperatures T = (main figure) and 



T = 50mK (inset), with Jmt(w) = corresponding to 
Rq — > oo. In the limit of low temperatures (fc^T /Hljqi <C 
1), the relaxation time becomes 



Ti(wqi) 
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(28) 



From the main plot of Fig. U several important points 
can be extracted. First, in the low frequency regime 
(woi "C fi) the RL (Drude) and RLC environments pro- 
duce essentially the same relaxation time Ti t jn,c(0) = 
Ti,rl{0) = Ti j0 ~ (Li/L) 2 RCo, because both systems 
are ohmic. Second, near resonance (woi ~ fi), T\,r.lc is 
substantially reduced because the qubit is resonantly cou- 
pled to its environment producing a distinct non-ohmic 
behavior. Third, for (<Joi > O), T± grows very rapidly 
in the RLC case. Notice that for wqi > V2Q, the RLC 
relaxation time Tx^rlc is always larger than Ti^rl. Fur- 
thermore, in the limit of <Joi S> max{f2,2r}, Ti^rlc 
grows with the fourth power of Wqi behaving as 



Ti 



RLC 



rp "01 



(29) 



while for woi 3> fl 2 /2T 1 T\ rl grows only with second 
power of (jJqi behaving as 



Ti 



RL 



4Ti. 



r 2 ^o 2 i 



(30) 



Thus, Ti t RLC is always much larger than Ti^rl for suf- 
ficiently large woi- For parameters in the experimental 
range (Fig. HJ, Ti^rlc is two orders of magnitude larger 
than Ti^rl, indicating a clear advantage of the RLC en- 
vironment shown in Fig Q] over the standard ohmic RL 
environment. Thermal effects are shown in the inset of 
Fig. 21 where T = 50mK is a characteristic experimental 
temperature 2 ^. Typical values of T\ at low frequencies 
vary from lCT 5 s at T = to l(T 6 s at T = 50mK, while 
the high frequency values remain essentially unchanged 
as thermal effects are not important for fiwoi 3> ksT. 

In the preceeding analysis we neglected the effect of 
the local environment by setting Yi nt (uj) — 0. As a re- 
sult, the low- frequency value of T\ is substantially larger 
than obtained in experimen t 12 ' 21 . By modeling the lo- 
cal environment with Rq = 5000 ohms and Lq = 0, we 
obtain the T\ versus wrji plot shown in Fig. [5] Notice 
that this value of Rq brings T\ to values close to 20ns 
at T = 0. The message to extract from Figs. [4] and [5] 
is that increasing Rq as much as possible and increasing 
the qubit frequency ujqi from 0. 10 to 2f2 at fixed low 
temperature can produce a large increase in T\. 

Having discussed the decoherence properties caused by 
environments with resonances for the case of supercon- 
ducting qubits, we discuss next how the environment 
shifts (renormalizes) the qubit frequency for the same 
systems. 
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FIG. 5: (Color-online) T\ (in nanoseconds) as a function of 
qubit frequency Woi- The solid (red) curves describe a phase 
qubit with RLC isolation network (Fig. |3]) with same parame- 
ters of Fig.[5]except that Ro = 5000 ohms. The dashed curves 
correspond to an RL isolation network with the same param- 
eters of the RLC network, except that C = 0. Main figure 
(T = 0), inset (T = 50mK) with ft = 141 x 10 9 rad/sec. 



IV. FREQUENCY RENORMALIZATION 

In this section, we discuss another effect of the environ- 
ment on qubit properties, in addition to dcphasing and 
relaxation. The environment also renormalizes (shifts) 
the qubit frequency cjoi by dressing the two-state system 
with environmental degrees of freedom. This is similar to 
the physics of the Lamb shift or the Franck Condon effect. 
In our case, the transition frequency E = uj$i is renormal- 
ized according to Er = E + 5E, where 5E = — Imi?ioio 
in terms of the Rcdfield rate tensor. The energy shift can 
be explicitly written as 

SE = —V \ duj , K ' [E coth(0w/2) - u] , (31) 

47T J E z - UJ A 

where V denotes the Cauchy principal value, and (3 = 
1/fcsT. Notice that 6E is analogous to the energy shift 
obtained in second order perturbation theory, which col- 
lects all processes in which a virtual boson is emitted and 
reabsorbed, such that no trace is left in the environment. 
The integral in Eq. (f3"Tj) can be calculated by extend- 
ing the integration to the complete real axis, closing the 
countour in the upper complex plane, and applying the 
residue theorem. 



A. Phase and flux qubits 

For phase and flux qubit discussed above, upon sum- 
mation over all residues of the relevant poles of the spec- 
tral density, we arrive to 2 ^ 
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K fl 3 E 
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G(r - ia) - RcG(iE) - n- 



E 



(32) 



where G(x) = + f3x/2ir) + ij}(f3x/2Tr) involves the 
digamma function ijj, and K = a($o/27r) 2 for phase 
qubit [sec Eq. (TJJ}], and K = a x [see Eq. flTDJl] for the 
flux qubit. Although, we use a compact notation involv- 
ing complex functions, the energy correction 5E is real. 
Notice that SE changes sign at u> ~ E, leading to an 
upward shift (renormalization) of E if most of the spec- 
tral weight of J(lo) is above E (corresponding to E < fl) 
whereas E shifts downward in the opposite case. Physi- 
cally, this corresponds to level repulsion between the spin 
and the oscillators in the environment. This result is con- 
sistent with usual second-order perturbation theory for 
energies. If the spectral weight J(oS) is concentrated at 
frequency u = f2, then the sign changes of SE happens 
at E ~ O, leading to a rather sharp structure in SE(Q,), 
and in Er(Q). 

In the limit of low temperatures (T — > 0), we can re- 
place the function G appearing in Eq. (|32[) by an appro- 
priate logarithm to find 



SE = 



K ifl 3 E 

27 r 



<T = ± 



e 2 - (an + ivy 



log 



r — iaQ 
iE 



(33) 

where £1 = £1 — T 2 /f2. In the underdamped limit il ^> T, 
one approximates the logarithm by log|£l/i?| — iair/2 
and rewrite the result as Er = E + SE, with SE = 
SEfi + SE ICS . The first term of SE contains a logarithmic 
contribution which resembles the scaling in the Ohmic 
case (with cutoff frequency fi), 



SE n = -J(E) log 

7T 



(34) 



where O w fi. This contribution changes sign from an 
upward shift at E > fi to a downward shift at E < fi as 
expected from the general arguments described above. 

The other contribution to SE takes into account the 
enormous spectral weight of the resonance, and can be 
written as 



5E ies = J{E) 



£ 2 -fi 2 

fir 



(35) 



A term of this kind persists even in the absence of damp- 
ing of the external oscillator. This term vanishes linearly 
with E at low energies, but undergoes the expected sign 
change near resonance E w fi, in which vicinity, a sub- 
stantial renormalization also occurs. 

As an illustration of the qualitative results discussed 
in this section, we show in Fig. [S] the frequency shift 
(renormalization) of the phase qubit with RLC isolation 
network described in Fig. [3] We make the identification 
E = ujqi and SE = Sujqi. Near resonance woi ~ fi, we 
find a frequency renormalization of about 2% which is 
due to the term 5E res . 

Next, we discuss briefly the frequency renormalization 
of charge qubits due to environmental effects. 
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FIG. 6: Renormalization of energy splitting for the phase 
qubit with RLC isolation network (Fig. [3]) for the parameters 
R = 50 ohms, Li = 3.9nH, L = 2.25pH, C = 2.22pF, and 
qubit parameters Co = 4.44pF, Rq = 5000 ohms and Lo = 0, 
T = 0, and Q = 141 x 10 9 rad/sec. 



B. Charge qubits 

For charge qubits with spectral density given by 
Eq. (|13[) . the energy renormalization 8E obtained in 
Eq. (Unj can be also calculated using complex integra- 
tion techniques. In the low temperature limit (T — > 0) 
there are two contributions to 



SE = SE n +SE res . 
The first one is a resonant contribution 

2 Clj A 2 [tt + 2 arctan(2fi 2 /r 2 )] 



SE r 



k 



where the function 



A 2 = e 1 - m 



and the second one is a non-resonant contribution 



SEa = 



2 cil r 2 » 2 io g [E 2 /(rl/4 + nj)] 



lc 



(A 



(36) 



(37) 



(38) 



(39) 



Just like the phase and flux qubit, the frequeny renormal- 
ization term due to the spectral weight of the resonance 
5E res is much larger than the logarithmic contribution 
SEq. While SE res undergoes a sign change from nega- 
tive to positive at E = y/Cl\ — (r 2 /2) 2 , SEq undergoes 
a sign change at yO| + (r2/2) 2 . 

Having concluded the discussion of the frequency 
renormalization (shift) of superconducting qubits due to 
coupling to environments with resonant frequency CI, we 
discuss next the behavior of superconducting qubits cou- 
pled to the same environments, when the qubit frequency 
loqi is close to CI, where essentially exact solutions are 
possible. 



NON-MARKOVIAN BEHAVIOR NEAR 
RESONANCE 



V. 



The Bloch-Rcdficld equations described in Eq. (|21j) 
capture the long time behavior of the density matrix, 
but can not describe the short time behavior of the sys- 
tem in particular near a resonance condition woi ~ CI, 
where the environmental spectral density J(tu ~ CI) is 
very large. In this case, only the environmental modes 
with cjfe = CI couple strongly to the two-level system, 
like a two-level atom interacting with an electromagnetic 
field cavity mode that has a finite lifetime. Thus, next we 
derive the time evolution of the state of the flux, phase 
and charge qubits when the qubit frequency is close to 
an environmental resonance. 

First, we restrict the Hamiltonian described in Eqs. §5§ 
and (IB]) only to boson modes with ujk ~ CI « cjoi- When 
dj^fiaiwoi, the Hamiltonian shown in Eqs. (0 and © 
can be solved in the rotating wave approximation us- 
ing the complete basis set of system-bath product states 

IV'o) = |0) s <8 |0) B ; |Vi) = |l)s ® |0) B ; |V*> = l«)s ® |*) B) 
where |0)s and |l)g are the states of the qubit and 
\k)s are the states of the bath. Notice that the states 
|l)s <S> |fc)B ar e absent in the basis set within the rotat- 
ing wave approximation and that the state of the total 
system at any time is 



4>(t) = CQ-ipa + ci(£)^i 



fc^O.l 



(40) 



with probability amplitudes cq, c\(t), and Ck(t). The 
amplitude Co is constant, while the amplitudes c±(t) and 
Cfc(t) are time dependent. Assuming that there are no 
excited bath modes at t = 0, we impose the initial con- 
dition Cfc(0) = 0, and use the normalization |</>(£)| 2 = 1 
to obtain the closed integro-differential equation 



ci(t) 



dttfit-tjaih), 



(41) 



where the kernel is the correlation function 

poo 

/(r) = / duj J(u>) exp [i(w i - w)r] 
Jo 

directly related to the spectral density J{u>). The reduced 
density matrix 



P[ ) ~ \ 4(t)c |c | 2 + E fe |c,(i)| 2 



(42) 



is subject to the condition Tr p(t) = 1, or more explicitly 
| Co | 2 + J2k \ c k(t)\ 2 — 1 — l c i(t)| 2 > which shows that the 
time dynamics of p(t) is fully determined by c\{t) for a 
specified value of cq. 



A. Phase and Flux qubits 

For phase and flux qubita with spectral densities given 
by Eq. (fT4"|) and Eq. (jTUJ) respectively, we can rewrite the 
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spectral density as 



Jres(w) = "777^ 2^ 



4ir 



<r=±l U) Z 



crfi + ir 



(43) 



where K = a($ /27r) 2 for phase qubit [see Eq. (fl"4"] . 
and if = «i [see Eq. (fTUj) ] for the flux qubit. This re- 
veals a resonance at lo — Q with linewidth T, imply- 
ing that any internal off-resonance contribution to J(u>), 
such as Jint{<^ = fi) in the case of the phase qubit cou- 
pled to an RLC environment, can be neglected for any 
non-zero value of the resistance Rq. In this case, the 
spectral density J{uj) dominated by the resonant contri- 
bution J(lu) w Ji SO (ui) « J Tes (uS). It is very important to 
emphasize, that evcnthough J Tes (uS) is rewritten explic- 
itly in terms of its poles in the complex plane, a simple 
inspection of Eq. (|4"3")l shows that J rcs (uj) is real. 

For this spectral density, we can now solve for c\{t) 
exactly and obtain the closed form 



c 1 (t) = £- 1 [ 



(s + r -iu 01 ) 2 + o 2 -r 2 



s [(s + r - iu ai ) 2 + n 2 - r 2 ] - Kn A ni/r 



where L 1 {F(s)} is the inverse Laplace transform of 
F{s), and k — K x ($ /2tt) 2 a l/(w iTi,o)- 

In Fig. [7J we plot the density matrix element pn = 
|ci(t)| 2 as a function of time for the dc-SQUID phase 
qubit described in Fig 03 The plot contains curves 
for three different values of resistance, assuming that 
the qubit is initially (time t = 0) in its excited state 
where pn(0) = 1. We consider the experimentally rel- 
evant weak dissipation limit of T <C u>oi ~ f2. Since 
r = 1/(2CR) the width of the resonance in the spectral 
density shown in Eq. (|43)) is smaller for larger values of 
R. Thus, for large R, the RLC environment transfers en- 
ergy resonantly back and forth to the qubit and induces 
Rabi-oscillations with an effective time dependent decay 
rate y(t) = ~2TZ{di{t) / Cl (t)} . 

These environmentally-induced Rabi oscillations are 
a clear signature of the non-Markovian behavior pro- 
duced by the RLC environment, and are completely ab- 
sent in the RL environment because the energy from the 
qubits is quickly dissipated without being temporarily 
stored. In the RL environment the decay in time of pu(t) 
has the characteristic non-oscillatory Markovian behav- 
ior. These environmentally-induced Rabi oscillations are 
generic features of circuits with resonances in the real 
part of the admittance. The frequency of the Rabi os- 
cillations Qrh = \J t^kQ 3 /2T is independent of the resis- 
tance since H,R a ~ ^l\J ~kL 2 C j L 2 Co, and has the value 
of fl Ra = 2nf Ra w 360 x 10 6 rad/sec in Fig. CO This 
effect is similar to the so-called circuit quantum electro- 
dynamics which has been of great experimental interest 
recently 2 ^— 

Next, we discuss briefly the non-Markovian dynamics 
for charged qubits. 



Pn 




200 400 600 800 1000 



FIG. 7: (Color-online) Population of the excited state of the 
dc-SQUID phase qubit (Fig. [3} as a function of time pn(i), 
with pn(t = 0) = 1 for R — 50 ohms (solid black curve), 
350 ohms (dotted blue curve), and R = 550 ohms (dashed 
red curve), and Li = 3.9nH, L = 2.25pH, C = 2.22pF, C = 
4.44pF, R a = oo and L = 0. 



B. Charge qubits 

For charge qubits, the spectral density in Eq. (fT3| can 
be rewritten in terms of its poles in the complex plane as 



e 2 n 2 

ile 



E 

cr=± 



lo — Q 2 + criT 2 /2 



(44) 



Again, we can solve for c\ (t) exactly, which is now given 
by 

d(t) = Cl (0)e- r ^/ 4 fcosh^ + ^smh^j (45) 



where A = - 8ne 2 fl 2 /lc. Notice that the fre- 

quency A can be complex, so it is convenient to rewrite 
it as 




where il c = Sire 2 jlc. In the underdamped case where 
^2 ^S> r 2 /2, provided that 2£l c /T 2 is not too small, A 
becomes purely imaginary: 



a w iy/n c n 2 . (47) 

Since cosh(i#) = cos(#) and smh(i6) = ism(0), where 
6 = (\/n c n 2 )t/2 in the present case, it becomes clear 
that an oscillatory decay of the state emerges in the un- 
derdamped regime. In this regime, the period of oscilla- 
tion is just t = 4ir/\/fl c n 2 . 

Having completed our discussion of environmentally- 
induced Rabi oscillations and its markedly 11011- 
markovian characteristic, we are ready to conclude. 

VI. CONCLUSIONS 

In conclusion, we analyzed decoherence effects in 
qubits coupled to environments containing resonances in 
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their spectral function, and we identified a crucial role 
played by the design of isolation circuits on decoherence 
properties. Furthermore, we found that the decoherence 
time of qubits can be two orders of magnitude larger 
than their typical low-frequency ohmic-rcgimc, provided 
that the frequency of the qubit is about two times larger 
than the resonance frequency of the environmental reso- 
nance (isolation circuit in the phase qubit case). We also 
studied the frequency renormalization (shift) of the qubit 
described by a two-level system due to dressing of the en- 
ergy levels by the environmental degrees of freedom. We 
found that the frequency shift changes sign across the res- 
onance frequency and is largest at resonance (about 2%). 
Lastly, we showed that when the qubit frequency is close 
to a resonance of the environment, the non-oscillatory 
Markovian decay of the excited state population of the 
qubit, gives in to an oscillatory non-Markovian decay, 
as the qubit and its environment self-generate Rabi os- 
cillations of characteristic time scales shorter than the 



decoherence time. In particular, we discussed as a con- 
crete example the decoherence properties of a dc-SQUID 
phase qubit coupled with an environmental RLC circuit 
possesing a resonance in its spectral function, where num- 
bers compatible with current experiments were used to 
estimate the environmental effects on decoherence prop- 
ertie a 21 ' 25 . 
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